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Dy 
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volume 
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Subscripts 
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TT 
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becomes  unstable 
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at  transition 

at  surface 
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INTRODUCTION 


The  object  of  this  investigation  is  to  develop  a  method  for 
predicting  the  location  of  transition  from  laminar  to  turbulent 
boundary-layer  flow  on  conical  re-entry  bodies.  Accurate  prediction 
of  the  transition  location  is  important  because  the  heat  flew  to 
a  surface  is  much  larger  when  the  boundary  layer  on  the  surface  is 
turbulent  than  when  it  is  laminar.  Because  of  the  high  heat-transfer 
rates  during  re-entry,  the  difference  between  turbulent  and  laminar 
heat- transfer  rates  is  so  large  that  the  variation  of  transition 
location  with  altitude  is  an  important  factor  in  the  choice  of 
surface  material  and  its  thickness  distribution  over  the  body. 

Material  that  is  too  thick  may  fail  because  of  thermal  strain. 

Material  that  is  too  thin  may  fail  by  burning  through.  Moreover, 
the  nose  shape  for  an  ablaLing  body  depends  on  the  variation  of  the 
transition  location  with  altitude  and  thus  so  do  the  aerodynamic 
characteristics  of  the  body. 

Attempts  to  predict  transition  location  began  about  forty  years 
ago  and  pertained  to  wings  and  bodies.  The  main  parameter  was 
usually  a  quantity  at  the  transition  point,  ordinarily  a  boundary 
layer-thickness  Reynolds  number.  The  assumption  that  transition 
occurs  at  a  constant  value  of  a  boundary-layer  Reynolds  number  some¬ 
times  gave  usable  results  for  a  more-or-less  fixed  pressure 
distribution.  Later,  Michel  (Ref.  (1))  developed  a  more  general 
and  seemingly  more  reliable  empirical  relation  between  Rgq  and  R£x 
at  transition  on  airfoils.  This  correlation  is  claimed  to  give 
sufficiently  accurate  calculated  profile  drag  coefficients  (Ref.  (2)). 
More  complicated  methods  nave  also  been  tried.  One  is  a  correlation 
of  (RS0t  -  Re0^)  with  the  average  value  of  the  velocity  profile 

shape  parameter  between  X>j  and  Xj_  (Ref.  (3)).  Another  (F.ef.  (4)) 
uses  a  family  of  similar  profiles  to  represent  the  velocity  profiles 
in  a  laminar  boundary  layer.  For  each  of  the  similar  profiles  the 
amplification  rate  of  a  small  oscillation  in  the  velocity  in  the 
laminar  boundary  layer  is  computed  as  a  function  of  the  Reynolds 
number.  The  amplification  rates  for  the  velocity  profiles  along 
the  boundary  layer  are  used  to  calculate  the  total  amplification  of 
a  small  disturbance  that  begins  at  the  stagnation  point  and  moves 
downstream.  The  total  amplification  up  to  the  measured  transition 
point  was  calculate^  for  a  number  of  cases.  Transition  usually 
occurred  for  an  amplification  of  between  about  e7  and  e12.  A  value 
of  about  e10  seemed  to  give  a  good  correlation  between  the  measured 
and  calculated  values  of  Rex  .  '-’his  method  is  complicated  and  seems 

not  to  be  used  by  Smith,  its  developer.  In  Smith's  latest  paper 
(Ref.  (2))  he  uses  Michel's  method.  Thus,  about  forty  years  after 
the  first  attempts,  the  preferred  method  for  the  prediction  of 
transition  in  incompressible  flow  seems  to  be  Michel's  method,  an 
empirical  relation  between  Re^  and  ReXT/  whose  limits  of  applicability 

are  indefinite.  Moreover,  M’chel's  relation  (Ref.  (2))  has  not  been 
shown  to  apply  when  there  is  flow  through  the  surface  or  when  the 
surface  is  rough. 
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For  compressible  flow  it  is  more  difficult  to  predict  the 
transition  location  than  it  is  fcr  incompressible  flow.  ^or 
incompressible  flow  over  smooth  solid  surfaces  only  the  pressure 
distribution  and  the  Reynolds  number  matter.  For  compressible  flow 
there  are  also  the  Mach  number  and  the  wall-to-stream  temperature 
ratio.  In  practical  problems,  flov  through  the  surface  and  roughness 
are  additional  parameters. 

As  already  mentioned,  the  ooject  of  the  present  investigation  is 
to  develop  a  method  for  predicting  the  transition  location  on  conical 
re-entry  bodies.  The  method  consists  of  attempting  to  find  variables 
that  allow  a  correlation  between  them  and  the  transition  location. 

Only  flight  and  ballistics-range  data  are  used  i  i  order  to  avoid  the 
transition-causing -disturbances  emitted  from  wind-tunne.l  wall 
turbulent  boundary  layers.  Although  some  sound  disturbances  can 
be  present  in  a  ballistics  range,  the  inference  from  Reference  5  is 
that  such  disturbances  probably  have  even  less  effect  on  transition 
than  those  in  a  wind  tunnel  whose  walls  are  covered  by  laminar  layers. 

ANALYSIS 

A  suggestion  for  choosing  correlation  parameters  is  sought  from 
the  velocity  fluctuation  kinetic  energy  equation  (Ref.  (C)).  An 
approximation  to  the  production  of  fluctuation  energy  is  taken  as 


-puTvT  and  an  approximation  to  its  dissipation  as  y (1-^— ) 

The  ratio  =•  between  production  and  a-  ;ipation  of  fluctuation  energy 
is  thus 


P 


D 


• — I — IT  3u 

-pu^vr  -T— 
.  aY 

,3u' . 2 

u(5y  > 


(1) 


P 

To  put  in  a  more  useable  form  the  terms  in  Equation  (1)  are 
written  in  non-dimensional  form  so  that  Equation  (1)  becomes 


(2) 
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Then  Equation  (2)  becomes 
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It  is  now  assumed  that  both 


*  u  ^ 


u '  v 


and 


to  the  maximum  value  of  ( — ; 

ue 


ue 


rD  H_. 
ue 


,0  I 


are  proportional 


(4} 


at  a  cross  section  of  the  boundary 
layer  and  thus  to  one  another.  Then  Equation  'O  becomes 


(v£)  (M£)  _f- 


(t) 


The  friction  coefficient  in  Equation  (5)  is  approximated  by  the 
reference  temperature  formula  for  a  laminar  boundary  layer  (Ref.  (7)), 


T 


f  _  k  ,u*.  ,p* 

"  R—  {~e]  {re) 
fc0 


The  quantity  v  in  Equation  (5)  is  evaluated  at  the  reference 
temperature.  Then  Equation  (5,  becomes 


(6) 
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(7) 


The  viscosity  is  assumed  to  be  given  with  sufficient  accuracy  by 


M£  _  To  111 
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Then  Equation  (7)  becomes 
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The  equation  for  the  rate  of  change  of  fluctuation  energy  is 
approximated  by 


peue 

or  with  Equation  (1) 

,2 


i  2 


bq*  ,3u  *  >  2  rP  ,  , 
peue  U  (3y  )  ID  -  11 


(10) 


(ID 


When  the  variation  of  ue  with  x  is  neglected,  Equation  (11)  can 
be  v bitten  as 


ua2 
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v7 
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The  ratio  (j^-)  in  Equation  (12)  is  replaced  by  (~)  and  Equation 
(9)  is  used  for  Then  Equation  (12)  becomes 
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It  is  assumed  that 


Te  Te  2 

is  negligible  when  compared  with  (^)  (^ )  ; 
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Equation  (14)  then  becomes 
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The  term 


3  — 

ue 

,3^/0 


is  assumed  to  be  proportional  to 


of  the  fluctuations  in  velocity  and  thus  to 
(15)  becomes 


the  magnitude 
Then  Equation 
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'•Tw;  VT*; 
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If  A  (^t)  has  to  reach  some  more-or-less  fixed  value  for  transition 
to  occur.  Equation  (17)  implies  that 


dx 
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For  the  purpose  of  this  rough  analysis  the  ratio  -r—  in  Equation 

dRe 

(18)  is  replaced  by  - — —  ;  this  neglects  the  variation  with  X  of 

Kea 

up  t 

the  local  Reynolds  number  per  foot,~.  Then  Equation  (18)  becomes 
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For  cases  in  which  the  integrand  in  Equation  (19)  varies  with  Rex  i 


roughly  the  same  manner,  the  quantity  [  (^)  (^)  ReeJ  should 

fall  in  a  more  or  less  narrow  band  when  plotted  against  ReXT*  That 


is,  approximately, 
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'eXn 


Tw  m  T#  2  _m 
F  (  (jj^)  (^)  ReQJi 


(20) 


The  value  2/3  for  m  in  Equation  (20)  is  assumed  to  be  precise 
enough  for  the  present  analysis.  Also,  the  expression  for 

T* 

(^■-)  in  Equation  (20)  is  taken  from  Reference  (7)  and  is 


T*  _  1  ,Tw  ,  ,  >  ,  v— 1  ,,  > 

Te  2  Te  +  15  +  *22r  2~  Me 


(21) 


For  the  present  correlation  attempt,  the  recovery  factor  r  in 
Equation  (21)  is  taken  as  .84  and  y  is  taken  as  1.4.  Equation  (21) 
then  becomes 


T*  _  1  ,Tw 
Te  “  2 


+  1)  +.037  Me2 


The  result  then  is 


Re*  =  F 


j‘S> 


'/3rl,Tw 
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2l4/3 


Re 
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(22) 


(23) 


The  data  used  to  test  the  supposition  (Eq.  (23))  are  those  in 
References  (3)  and  (9) .  The  data  from  Reference  (8)  are  for  no 
ablation.  There  are  three  tables  of  data  in  Reference  (8) .  The 

R~  /  Rea  i  and  Mem  which  were  used  in 
XT  Tw  eT  1 

of  (=— )  _  needed  in  Equation  (23)  were 
Tw  ^  ^  ^ 

obtained  from  the  values  of  (=— )  given  in  the  tables  of  Reference  (8) 
by  the  expression  T 


tables  give  the  values  of 
Equation  (23)  .  The  values 


H  -  gfU+/PF  Ke2), 


(24) 


or  witn  /Pr  =  .84  and  y=1.4, 


(1+ .  168  Me2) 


(25) 


The  first  table  in  Reference  (8)  is  headed,  "Temperature  ratio 
investigations."  This  table  contains  data  for  sharp  5-degree  half¬ 
angle  cones  tested  through  a  f roe-stream  Mach  range  from  2.95  to 
8.29.  The  second  table  is  headed,  "Mach  number  and  body  geometry 
investigations."  This  table  contains  data  for  3t  5y  6.37  and  9-degree 
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half-angle  cones  with  riose-to-base  radius  ratios  from  .012  to  .10. 

The  f ree-stream  Mach  numbers  varied  from  6.7  to  12.8.  The  third 
table  is  headed,  "Drag  Method."  This  table  contains  data  for  6.3- 
and  9-degree  half-angle  cones  with  nose-to-base  radius  ratios  from 
.017  to  .12.  The  f ree-stream  Mach  numbers  ranged  from  9  to  15.  The 
Tw 

range  in  covered  by  the  complete  set  of  tests  listed  in  the  three 
tables  was  from  .026  to  .265. 

The  parameter  {  <~)  2/3  lj  +  1)  + .  037Me2  J  4/3  Re^  )  ,  called 
G„, ,  is  plotted  against  Rev  in  Fiqure  1(a);  each  set  of  data  has 

I  Arp 

its  own  symbol  (see  Table  1) .  The  data  lie  in  a  fairly  narrow  band. 

The  variation  in  Re.<T  for  the  same  GT  is  in  general  no  larger  than 

the  variation  in  RSx^  Reference  (8)  for  the  same  flight  conditions. 

The  data  shown  in  Figure  1(a)  were  fitted  by  a  least-squares 

equation  of  the  type 

P.fixT  =  a  0/  (26) 

with  the  result 

=  2.259  x  JO4  G  *7764  (27) 

*-XT 

The  plot  of  this  equation  is  shown  in  Figure  1(b).  The  Equation 
(27)  it.  plotted  on  a  separate  figure  from  Figure  1(a)  so  that  the 
distribution  of  data  in  Figure  1(a)  can  be  observed  without  the 
intrusion  of  the  effect  of  a  curve  through  the  data.  The  dashed 
portion  of  the  curve  in  Figure  1 (b)  is  an  extrapolation  beyond  the 
range  of  the  data.  The  Equation  (27)  seems  to  fit  the  data  with 
acceptable  closeness.  A  better  representation  can  probably  be 
obtained  by  using  a  more  complicated  expression  than  Equation  (26)  or 
even  by  drawing  in  a  curve  by  "eye." 

In  Figure  2(a)  the  40  points  of  "high-quality1'  flight  data 
from  Reference  (9)  are  shown.  The  identification  code  for  each 
point  in  Figure  2(a)  is  given  in  Table  II.  The  data  cover  a  range 
of  free-stream  Mach  number  at  transition  from  3.17  to  23.70.  The 
data  include  a  range  of  cone  angle,  nose  radius,  ana  various 
combinations  of  materials  for  the  nose  and  frustum  of  the  cones. 
Specific  information  is  given  in  a  classified  appendix  to  Reference 
(9)  . 

When  the  points  G^.,  Re  were  first  plotted,  those  for  large 
•  ^ 

values  of  mf,  the  ablation_on  the  frustum  of  the  cones,  fell  too 

far  from  those  with  small  for  the  present  method  of  correlation 

•  • 

to  be  satisfactory.  The  quantity  is  called  m  in  Table  7 
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of  Reference  (9) .  It  was  found  that  a  satisfactory  correlation  for 
all  the  40  points  of  Reference  (9) ,  both  those  for  high  and  low 
ablation,  was  obtained  when  Gt  was  replaced  by  G*T  where 


G*<r  ~ 


and 


g*t  -  —= 


mf 


-3 


for  £  10 


x  10  3  for  ~  10  3, 


The  correlation  shown  in  Figure  2(a)  is  G*p  against  RGxT.  The  values 

of  iii£  are  those  given  in  Table  7  of  Reference  (9)  .  The  values  of 

Tw  Tw 

(=— )  needed  in  Equation  (23)  wore  obtained  from  the  values  of  (=— ) 
rp  ^  T 

given  in  Table  7  of  Reference  (9)  by  the  relation 


Tw  _  Tw 
Te  T«> 


1  ♦  ^  Me2 

1  +  A)—  Moo 2 


(28) 


or,  for  y  *  1.4, 


Tw  .  Tw  5  -i  Me2 
Te  T«,  5  +  Moo2 


(29) 


A  least-squares  expression  tnat  fits  the  data  of  Figure  2(a) 
in  a  satisfactory  manner  is 

Rg  =  3.120  x  104  G*"*72?<  (30) 

A  plot  of  this  expression  is  shown  in  Figure  2(b). 

Figure  3(a)  shows  all  the  data,  both  ballistics  range  and 
flight.  A  least-squares  expression  that  gives  a  satisfactory  fit  for 
ali  the  data  in  Figure  3(a)  is 

Re  =  2.932  x  104  G*T  *7370  (31) 

A  plot  of  this  expression  is  shown  in  Figure  3(b). 
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Finally,  summing  up,  the  expression  which  correlates  all  the 
data  considered  is, 


ruxt 

®  2 

.932 

X 

104 

G*t‘ 

7370 

(32) 

where 

G*  = 

G 

for 

mf  i 

10-3  ) 

and 

A 

(33) 

G*  - 

G 

x  10 

-3 

for 

"f  - 

10"3 

*f 

and 

4 

G  = 

(— ) 
lTe; 

2/3 

l7 

,Tw 

Te 

+  i) 

+  .037 

Me2]7  EwC 

(34) 

The  data  correlated  by  Equation  (32)  cover  a  range  of  MeT  from 

Tw 

2.8  to  14.7,  a  range  of  (=—)  from  .074  to  1.43,  and  a  range  of 

—  1  e  T 

•  1  -3 

ablation  parameter  from  0  to  17.7  x  10 


DISCUSSION 

It  is  obvious  from  the  derivation  of  the  relation 


Rexf  =  ^ 


(35) 


that  many  approximations  have  been  made.  Consequently,  although 
a  satisfactory  correlation  has  been  obtained  for  the  88  points  in 
Figure  3(a),  there  is  no  guarantee  that  experiments  in  which  some 
conditions  differ  sufficiently  from  those  for  these  88  points  will 
yield  values  of  G*T  and  Rex  will  conform  to  the  present 

*  Tw 

correlation.  One  point  that  must  be  stressed  is  that  Rex»  Me» 

and  Reg  must  be  calculated  by  the  same  procedures  used  in  References 
(8)  and  (9).  The  calculation  procedures,  especially  for  blunt  cones 
and  for  cones  whose  shapes  are  changing  because  of  ablation,  are 
part  of  the  correlation  procedure.  If  other  calculation  methods  are 
used,  the  user  must  be  prepared  to  obtain  new  correlations,  hopefully 
involving  the  same  variables  used  in  the  present  investigation. 


In  Michel's  method  (Refs.  (1)  and  (2))  the  computed  curve  of 
Reo  versus  Rex  intersects  the  correlation  curve  of  ReOT  versus  Rexx 

at  a  small  angle.  This  often  makes  it  difficult  to  obtain  a  definite 
value  for  Keo_  or  Re  A  trial  of  the  present  procedure  was  made 

for  shot  number  5885  and  shot  number  5874  (Table  I),  cc  ics  for  which 
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Te  and  Me  varied  with  X.  For  each  of  these  two  cones  the  intersection 
angle  between  the  curve  of  G  against  Rex  and  the  correlation  curve  of 
Figure  3  was  large  enough  so  that  Rex  could  be  read  without  difficulty. 

This  is  probably  a  general  result  for  the  present  method.  Thus,  for 
Tw 

and  Me  independent  of  X,  Equation  (32)  becomes 

R**T  *  C°  %T'7370  (36> 

For  a  sharp  cone,  however, 


Rex  «  C1  Reg 2  (37) 

The  difference  in  the  value  of  the  exponent  of  Reg  in  Equations 
(36)  and  (37)  probably  implies  an  angle  of  intersection  large  enough 
to  eliminate  difficulty  from  this  source. 

The  roughness  of  the  surface  does  not  appear  as  a  parameter  in 
the  expression  for  G*  even  though  transition  is  known  to  be  affected 
by  roughness.  Eecause  a  correlation  of  the  88  points  is  obtained 
without  a  roughness  parameter  it  probably  means  that  the  roughnesses 
of  the  various  materials  were  small  enough  not  to  exert  a  noticeable 
effect. 

p  v 

^  w  w 

Like  the  roughness,  the  blowing  or  ablation  parameter  — does 

PeJe 

not  appear  explicitly  in  G*.  The  calculated  Reynolds  number  Reg* 

p  v  p  v 

however,  depends  on - and  so  G*  does  depend  on  -  . 

peue  Peue 

Equation  (20)  was  obtained  from  Equation  (19)  by  limiting  the 
cases  to  those  for  which  the  integrand  in  Equation  (19)  has  roughly 
the  same  type  of  variation  with  Rev*  All  of  the  data  in  References 
(8)  and  (9)  seem  to  be  for  transition  on  the  cone  frustum.  Consequently, 
if  transition  occurs  on  a  blunt  nose  the  variation  of  the  integrand 
in  Equation  (19),  which  is  G"l,  with  Rex  may  differ  sufficiently  from 
its  variation  on  the  frustum  for  a  pointG*T,  ReXT  to  fall  too  far 

from  the  band  of  points  in  Figure  3.  Experiments  in  which  transition 
occurs  naturally  on  smooth  blunt  noses  are  therefore  desirable  to 
test  the  correlation  method's  ability  to  predict  transition  on  blunt 
nosec.  Some  modification  of  the  correlation  may  be  necessary. 

CONCLUSION 

An  empirical  equation, 

Re  =  2.932  x  104  G*  ‘73/0, 

^xt  1 

has  been  developed  which  predicts  the  transition  location  fairly  well 
for  48  ballistics  range  tests  for  cones  without  ablation  and  for 
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40  flight  tests  for  cones  with  ablation  ranging  from  zero  to  large. 

The  experimental  data  cover  a  range  of  MeT  from  2.8  to  14.7,  a 

Tw 

range  of  (=—)  from  .074  to  1.43,  and  a  range  of  ablation  parameter 

_  T 

mf  from  0  to  17.7  x  10-^.  Nose  tips  range  from  sharp  to  fairly  blunt. 
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TABLE  I  -  IDENTIFICATION  OF  BALLISTICS  RANGE  DATA  POINTS 


Temperature  Ratio  Investigation  X 


Shot  No. 

ReXT*  10'6 

Shot  No. 

Rex„*  10'6 
* 

4046 

?  "7 

5041 

5.25 

4047 

V.  • 

5042 

4.40 

4048 

...  " 

5043 

9.20 

4051 

2.35 

5049 

7.67 

4056 

6.68 

5050 

5.36 

4057 

6.04 

5467 

4.98 

4058 

5.81  j 

546S 

4.42 

4059 

3.96 

5  *69 

4.91 

4060 

3.13  ! 

5470 

4.69 

5036 

5.46  ! 

5462 

9.32 

5037 

5.25 

5474 

11.56 

5039 

9.23 

5475 

7.29 

5040 

8.11 

j 

! _ 

_ 

Mach  Number  and  Geometry  Investigation  + 

5697 

8.7 

5840 

1.3 

5699 

14,2 

5343 

11.1 

5707 

11.4 

5857 

12.6 

5708 

6.0 

5874 

0.5 

5776 

0.6 

5881 

5.3 

5779 

7.0 

5885 

15.2 

5802 

12.1 

5895 

1.9 

5836 

1.1 

Drag 

Method  Q 

D1 

6.8 

D5 

2.4 

D2 

6.8 

D6 

1.0 

D3 

16.3 

D7 

3.4 

D4 

5.6 

D8 

3.3 

Reference 

(8)  for  more  information  for 

each  data  point 

Data  points  that  do  not  appear  in  Figures  1  or  3  have  values  of 

gt<  100. 


NOLTR  73-127 


TABLE  Ii  -  IDENTIFICATION  OF  FLIGHT  DATA  POINTS 


Body 

Symbol 

HeXT  ^  10 

Symbol 

ReXT  X  10“6 

0 

49.4 

C05 

D 

12.3 

0 

29.1 

C06 

O 

4.35 

B03 

0 

21.3 

SJ1 

0 

51.8 

B04 

0 

18.8 

S02 

0 

86.5 

B05 

0 

15.6 

S03 

0 

47.4 

B06 

0 

7.17 

S04 

0 

34.6 

B07 

0 

30.3 

SO  5 

'7 

4.78 

B08 

0 

1.45 

S07 

9.  50 

BIO 

0 

0.81 

S  08 

0 

17.0 

B13 

0 

2.53 

S09 

0 

1.18 

B15 

0 

2.36 

S10 

0 

29.1 

BO  9 

n 

13.7 

T01 

0 

11.1 

Bll 

o 

0.40 

T02 

0 

41.1 

B12 

O 

0.  33 

T03 

0 

16.5 

B14 

O 

0.44 

T05 

0 

14.6 

B16 

O 

0.47 

T06 

0 

16.1 

C01 

V 

5.3 

T07 

0 

15.8 

C02 

21.4 

T08 

0 

15.3 

C03 

'll 

9.18 

T04 

24.1 

C04 

6.32 

T10 

_ 

_ 

12.5 

_ - 

See  Reference  (9)  for  more  information  for  each  data  point. 


Data  points  that  do  not  appear  in  Figures  2  or  3  have  values  of 

g*t  <100. 
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FIG.  3  (a)  VARIATION  OF  TRANSITION  PARAMETER  G-,  WITH  TRANSITION  REYNOLDS  NUMBER 
R«,  FOR  BALLISTICS  RANGE  AND  FLIGHT  DATA,  EXPERIMENTAI  DATA. 

(See  Tables  I  and  II) 


FIG.  3  (b)  VARIATION  OF  TRANSITION  PARAMETER  G*y  WITH  TRANSITION  REYNOLDS  NUMBER 
R  FOR  BALLISTICS  RANGE  AND  FLIGHT  DATA,  LEAST  SQUARES  EQUATION  FOR 
e*T 

EXPERIMENTAL  DATA. 


